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Abstract 



On , We derive contributions to the trace formula for the spectral density accounting for the role of 

— ..., I diffractive orbits in two-dimensional polygonal billiards. In polygons, diffraction typically occurs at 

^ ' the boundary of a family of trajectories. In this case the first diffractive correction to the contribution 

i-rt ■ of the family to the periodic orbit expansion is of order of the one of an isolated orbit, and gives the first 

J^ ! yh correction to the leading semi-classical term. For treating these corrections Keller's geometrical 

j^ ' theory of diffraction is inadequate and we develop an alternative approximation based on Kirchhoff 's 

1-^ ■ theory. Numerical checks show that our procedure allows to reduce the typical semi-classical error 

. . by about two orders of magnitude. The method permits to treat the related problem of flux-line 



diffraction with the same degree of accuracy. 
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1 Introduction 

Two dimensional billiards play a central role in the domain of quantum chaos because of the simpli- 
city of their classical dynamics and of the relatively easy determination of their quantum spectrum. 
During the last 20 years they have been used as model systems for testing semiclassical trace formulae 
(following Gutzwiller Q and Balian and Bloch Q) and random matrix theory (see e.g. Q). 

Amongst these systems, plane polygonal billiards have been subject of a long lasting interest (see 
e.g. the review Q) : they have zero metric [|[ and topological |^ entropy, but their dynamical 
properties range from integrable to possibly ergodic and mixing Q, passing by the interesting group 
of pseudo-integrable systems [^. Level correlation of integrable polygonal billiards display interes- 
ting properties g, not to speak of the case of pseudo-integrable billiards whose level statistics are 



intriguingly related to those of the Anderson model at the metal- insulator point [10, 11|. 

The present work is devoted to the detailed study of the trace formula in polygonal billiards. 
Though the general method of deriving the trace formula is well known ||l[ , its application to polygonal 
plane billiards is not straightforward. The main difficulty is the existence of important corrections due 
to the diffraction on the corners of the billiard. This type of correction was treated in Refs. p^-p!^] 
in the framework of Keller's geometrical theory of diffraction |15]. This amounts to introduce in the 
trace formula new diffractive orbits that obey the laws of classical mechanics everywhere except on 
singularities of the potential where they are diffracted non-classically. The result of the approach 
of Refs. [0-14 1 diverges when a diffractive orbit is close to be allowed by classical mechanics ; this 



deficiency was remedied in some special cases in Refs. [16, 17[. Ref. [jT^ studies corner diffraction in 
two dimensional billiards (not exclusively polygons). It gives uniform formulae but is limited to single 



diffraction. Ref. [jlgl treats diffraction by a circular disk inside a billiard. It considers up to doubly 
diffractive orbits, but does not provide a uniform approximation. In the present paper we extend these 
approaches and construct improvements to the geometrical theory of diffraction in polygonal billiards. 

This type of corrections is made necessary in polygons because in these systems the spatial ex- 
tension of a family of orbits is often stopped by a singularity of the frontier of the billiard ; as a 
result the generic situation is that a diffractive orbit appears on the boundary of each family. This 
trajectory is on the verge of being allowed by classical mechanics and thus cannot be included in the 
trace formula in the framework of the geometrical theory of diffraction. Hence in the following we 
devote a special care to the treatment of diffractive periodic orbits lying on the boundary of a family 
and of its repetitions. We give explicit formulae for the corrections to the leading semi-classical term 
for the n*'' iterate of a family of periodic orbits. 

We find in polygonal billiards a very rich variety of diffractive orbits. Their contributions give 
V^ corrections to the leading semi-classical term in the trace formula and allow to compute the level 
density with great precision. Numerical checks show that the typical semiclassical error is reduced by 
one or two orders of magnitude. 

The paper is organized as follows. In Section |2| we briefly present Keller's geometrical theory 
of diffraction and propose an alternative approximation based on Kirchhoff's theory that is valid 
near the "optical boundary" (the separation between allowed and forbidden classical trajectories, in 
optics this occurs on the line separating light and shadow). The simplicity of the method permits a 
straightforward generalization to the case of diffraction by a flux line. The approximation established 
in Section |2| is used to treat a large number of different types of diffractive periodic orbits. We first 
consider corner diffraction. The contribution of a diffractive orbit on the boundary of a periodic orbit 



family is calculated in Section ^. This is a typical situation for pseudo-integrable billiards. Special 
attention is given to the diffractive partner of the n-fold repetition of a primitive periodic orbit. Section 
^ is devoted to the study of diffractive orbits that are simultaneously on the boundary of a family 
and on the frontier of the billiard. Another type of diffractive orbits that belong to the boundaries 
of two different families of periodic orbits is discussed in Section ^. Besides diffractive orbits lying 
exactly on an optical boundaries, there exist orbits that are so close to an optical boundary that 
the geometrical theory of diffraction cannot be applied. Such orbits are studied in Sections ^ and 
^. All these special cases are necessary for a careful description of the quantum density of states in 
pseudo-integrable billiards. In Sec. ^ we illustrate the flexibility of our approach by presenting results 
for flux line diffraction in a rectangular billiard. In this case, solving the question of diffraction on 
the optical boundary amounts to treat the non-trivial problem of (multiple) forward Aharonov-Bohm 
scattering. Finally we present our conclusions in Section 0. Some technical points are given in the 



Appendices. In [Appendix A| a concise discussion of improvements to Keller's theory of diffraction is 



given. In Appendix B we discuss the computation of certain trace integrals. In Appendices C and D 



we derive analytically explicit expressions for important n-dimensional integrals. 

2 Diffractive Green function 

In this Section we first present Keller's geometrical theory of diffraction (putting the emphasis on 
corner diffraction) and then propose an alternative approximation valid near the optical boundary 
(when Keller's approach fails) for corner and flux-line diffraction. 

2.1 Geometrical theory of diffraction 

One considers the different approximate contributions to the Green function G{f, f , E) for two points 
r and f in a polygonal billiard. The first is the semi-classical contribution which is a sum over all 
possible classical trajectories going from f to r. It is of the form : 

G„(r. f, E) = Y. '°"'''"^^'''^'» , (1) 

classical 

where L is the length of the classical path going from f to r and u is twice the number of specular 
reflections along that path (we consider Dirichlet boundary conditions). We use units such that the 
energy is related to the wave- vector through E = k'^. 

There are other contributions to G that correspond to diffractive orbits experiencing specular 
reflections on the frontier of the billiard and also non-classical bounces on the diffractive corner. In 
the framework of Keller's geometrical theory of diffraction (see e.g. |jl5|) a such orbit with a single 
diffractive bounce contributes to the Green function with a term 



Gid(f, r ', E) = Go(fo, r ', E) V{e, 9') Go{r, fo, E) , (2) 



where rg is the position of the diffractive apex and T> is a diffraction coefficient depending on the 
interior angle 7 of the polygon at ro and on the incoming (outgoing) angle 6' (9) of the diffractive 
trajectory with the boundary. The explicit expression of V for corner diffraction reads ||I^ : 

Vie, 9')= J2 T^^AO^O') where P,,,(^,0') = ^ zfr^)^ ' ^^^ 

tan (^^^ 

where A^ = j/ir and ^o- = 0' — o"0, (with 9 and 9' in [0,7]). Although formula (Q) is attached to the 
name of Keller, the idea of treating diffraction as arising from a kind of reflection on the edge has a 



long history which goes back to Young (see Chap. 44 of |18| and Chap. 8.9 of [^). We also note here 
that the repercussion of diffractive periodic orbits on the spectrum seems to have been first worked 
out within the geometrical theory of diffraction by Durso in 1988 pO]. 

Formula (|2|) can be generalized to treat multiple diffraction. One has then several diffraction 
coefficients 2?i,I?2) ■ ■ ■ , one for each diffractive bounce, and between each diffractive bounce a semi- 
classical propagation described by a Green function of type (||). When diffractive trajectories are taken 
into account in the trace formula, one is lead to consider diffractive periodic orbits which contributions 



to the level density p{E) are of the form [12-14] : 




piE) ^lAiU ^=^ } cos(fcLd - i^d7r/2 - 3n7r/4) . (4) 



In (0) and in many instances below, when writing explicitly the contribution of a periodic orbit 
(classical or diffractive) to the level density, we put an arrow in direction of p{E) for indicating that 
this is one contribution amongst many others. In the above expression ii, ...,in are the lengths along 
the orbit between two diffractive reflections. £1 + ... + in = L^ is the total length of the diffractive 
periodic orbit, v^ is the Maslov index of the diffractive orbit, i.e. twice the number of specular 
reflections. Repetitions of a primitive diffractive orbit appear as a special case of (§) ; in this case 
however, in the first factor Ld/(27r/c) of the r.h.s. of (^), L^ should be understood as the primitive 
length of the orbit. 

We recall that in a polygon, the contribution of an isolated periodic orbit to p{E) is of the form 
— £/(47r/c) cos(A;£) (for a primitive orbit of length €). Thus Eq. (^) shows that the contribution of a 
typical diffractive periodic orbit with n diffractive bounces is of order 0{k^"''^) compared to that of 
an isolated periodic orbit. We will study below special configurations where this is not the case and 
where diffractive orbits have the same weight as isolated periodic ones. We first need to discuss the 
range of validity of the geometrical theory of diffraction and to define approximations alternative to 
(I)- 

2.2 In vicinity of an optical boundary 

The approximation defined by Eqs. (|2|,y) fails when the diffractive bounce at rg is "almost allowed" 
by classical mechanics. In that limit the trajectory lies on what is called an "optical boundary" in 



the literature, and the coefficient T> diverges. This failure of Keller's approximation can be intuitively 
understood by noting that Eq. (|2|) gives a contribution to the Green function of order 0{k~^) whereas 
in the limit that the diffractive orbit becomes allowed by classical mechanics it has to contribute at 
order ©(fc^^'^) as any classical trajectory. Hence in this limit Eq. (|2|) cannot hold. We study below 
a triangle with a diffractive corner of opening angle 7 = Stt/S. In this case, one can easily check 
geometrically that the diffractive orbit coincides with a classical trajectory if the angles 6 and 6' lie 
on one of the lines of [0, 7]^ shown on Fig. |l|. This can be also checked algebraically from formula (^) : 
each of the four lines of Fig. 1^ correspond to divergence of one of the coefficients 'Da,ri- 
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Figure 1: Solid lines : location of the angles 6 and 9' for which expression (ph of V diverges when 7 — Stt/S. 
Near each line are indicated the values a and 77 of the divergent Po-.r;- Dashed lines : location of the angles 9 
and 9' for which expression (||) of V is zero. These correspond to 9 (or 9') = or 7. 



For corner diffraction, after the work of Pauli |21|, several uniform approximations have been 
derived which correct the drawbacks of Eq. (^). We recall one of these in Appendix A.l| . In this 
paper we use a simple approximation to the exact formula valid only near the optical boundary. Let's 
consider that the trajectory lies near the optical boundary defined by one of the four couples (cr, r]) ; 
then our approximation for the total Green function (semi-classical plus diffractive) reads : 



r+co 

Gi(r, r', E) = -2 / ds Go(,s, f', E) n^ ■ V,-Go(f, s,E) 
Jo 

+ Go(fo, r ', E) Vregie, 6') Go(f, fo, E) , 



(5) 



where the locus of points s is an arbitrary half- line separating r and r and issued from fo (at s = 0); 
fis is a vector normal to the s axis and oriented towards f, see Fig. |2|a. P^eg is the non-divergent 
part of the diffraction coefficient (i.e. the sum of all the 2?o-,?7's but the divergent one). The diffractive 
Green function (analogous to (^)) is defined from Eq. (^) by the difference : Gid = Gi — Gq. 



Eq. (|^ is a simple Kirchhoff approximation to the Green function (with Keher-type corrections) 
which is to be used within the semiclassical approximation (this is ihustrated at length in the following 



sections). We show in Appendix A how it can be derived starting from a more elaborate approach. 
Eq. (g) is exact in the limit that the classical path from r ' to r lies on an optical boundary. It is 
designed to remedy the divergence of the geometrical theory of diffraction and it is not a uniform 
approximation to the Green function : far from the optical boundary characterized by a and rj, it 
yields a result basically of the form (g) without however the correct form of the coefficient ^o-,»? (this 
is in accordance with the known aspects of Kirchhoff approximation, see e.g. [|l8l). This should not be 
considered as a limitation of the approach : we show below (Sees. |6| and j^) that it is a simple matter to 
"uniformize" the result derived from (^) when necessary. Compared to the uniform expression (|Alj), 
Eq. (|^ has the important advantage of being easily extended to treat multiple diffraction near the 
optical boundary (Eq. @)- In the following sections approximation (^) will allow us to incorporate 
non-standard diffractive contributions in the trace formula. 



(a) 



b ) 





Figure 2: Graphical representation of the notations of Eqs. (H) and (o). Part (a) refers to corner diffraction 
and part (b) to flux line diffraction, tq is the diffractive point. The dashed line issued from vq is the optical 
boundary on which the geometrical theory of diffraction fails. In part (a) the integration along s is stopped on 
the apex at tq- This is not the case in part (b). There however, near the optical boundary \S(j)\ <ti n and the 
integrand of (j^ contributes with a phase that is approximatively —an if s > and an if s < 0. 

Formula (|5|) can be extended to treat the case of diffraction by a flux line. This problem bears 
important similarities with corner diffraction. In some respect it can be considered as simpler, because 
for an initial point r , the diffractive point (the Aharonov-Bohm flux line) is associated with a single 
diffractive boundary : the forward direction. This is the reason why formula @ below - which is the 
analogous of Eq. (||) - comprises only a Kirchhoff contribution and no Keller-like correction. 

We consider a particle of charge q and a flux line located on point fg such that the magnetic field 
is -B = $5(r — fo )z. The only relevant parameter is the ratio a of the flux $ with the quantum of flux 



6 



Q. = q(t/{2'Khc) and one can restrict oneself to < a < 1. The KirchhofF approximation for the total 



Green function is (see the derivation in Appendix A. 2 ) 



Gi(f, r ', E) = -2 f ds Go{s, r ', E) n, • VgGoif, s, E) e^«^<^(^) , (6) 

J —oo 

where the locus of points s is an arbitrary line separating f ' and r and going through vq (at s = 0) , 
and A(/>(s) is the angle covered by the path going from f to s and then to r. If Scj) is the angle 
between r — tq and rQ — f (i.e. the departure from the optical boundary) then A(/)(s) = 5</> — tt sgn (s). 
Of course, in this procedure, the orientation of the axis (ro,s) is not arbitrary. Our choice of A(/)(s) 
corresponds to an orientation such as presented on Fig. gb. 

Eq. (P) has a simple physical interpretation : the phase accumulated by a trajectory depends on 
the sense of rotation of the circuit around the flux line. As Eq. (g) it is only valid near the optical 
boundary, but it is easily generalized to multiple diffraction in the forward direction : i.e. it allows to 
treat the problem of multiple forward Aharonov-Bohm diffusion. We illustrate this property in Sec. ^. 

3 A diffractive orbit on the frontier of a family 

A typical occurrence of diffractive orbits is at the boundary of a family of trajectories. The width 
of a beam of classical orbits is limited by a singularity of the frontier of the billiard. Such a case is 
illustrated by the example shown in Fig. y. Note that this is not the only possible type of boundary 
of a family. It may happen that the family stops on a non-diffractive corner (a corner with opening 
angle of the type vr/n with n G N). This is the case for one of the boundaries of the family displayed 
in Fig. y. The frontier of the family may also be one of the frontiers of the billiard, this is illustrated 
on Fig. 1^. We will also study below (Section H) as mixed case where the boundary of the family only 
partly coincides with the frontier of the billiard. 

On Figure ^ we have represented the family by one of its member (upper left triangle). The 
boundary of the family is shown in the lower left triangle. Also, to the right of the plot, instead 
of representing the orbit by a series of segments bouncing off the frontier of the billiard, we have 
represented it by a unique straight segment where the reflection on each edge is replaced by continuing 
the path into a reflection of the enclosure. This procedure is called "unfolding the trajectory". 

The diffractive orbit on the boundary of a family appears as a correction to the contribution of 
the family and of its repetitions. Its contribution to the trace formula cannot be evaluated from the 
geometrical theory of diffraction, because its coefficient T> is infinite. However, since the diffractive 
orbit is exactly on an optical boundary, it can be described by using Eq. (^). The contribution 
of an orbit to the level density is evaluated in the framework of a semiclassical periodic expansion 
(see e.g. Refs. H, ^) : the level density being related to the Green function through the trace 
p{E) = — (l/vr) Im J (Pr G{f , f , E + iO^) , the Green function is approximated in vicinity of a periodic 
orbit (here it will be done by using Eq. (g) but other approximations will be used below) and the 
trace integral is evaluated within a saddle phase approximation near the saddle corresponding to the 
periodic orbit considered. 

• The diffractive correction to the first iterate of a family is very simple : the trace of the first 
diffractive contribution to the Kirchhoff Green function in (Eh is zero and only the term with a regular 




Figure 3: Representation of a periodic orbit which is part of a family in the triangle (7r/2, 7r/8,37r/8) (upper 
left triangle) and of the diffractive orbit which is on the boundary of the family (lower left triangle). The plot 
to the right represents the family after unfolding. The area occupied by the family is shaded and the diffractive 
point on its boundary is marked with a black dot (as in the lower left triangle). The diffractive orbit in the 
lower left triangle appears as the right boundary of the unfolded trajectory. 

Keller-type diffraction appears. Hence, if the family has a length d. and occupies on the billiard an 
area A, then its total contribution (semi-classical plus diffractive) to p{E) is simply : 

p{E) ^ ^^i= cos{M - 7r/4) + ^/r^f^ cos(A:£ - v,tt/2 - 3vr/4) . (7) 

The first part of the r.h.s. of (0) is the usual contribution of a family of periodic orbits in two 
dimensions. The second part is of the form (H) : it comes from the regular Keller contribution in (0). 

In order to test the validity of our approach, we have compared our analytical results with the 
spectrum determined numerically in a triangular billiard with angles (7r/2,7r/8, Svr/S). Note that 
this is not a generic polygonal billiard : its classical mechanics is pseudo-integrable and furthermore. 



it belongs to the set of "Veech billiards" |22]. We choose these systems because they simplify the 
geometrical computations : Veech billiards have the amusing property that there exists only a finite 
number of possible areas occupied by a family of periodic orbits. In the triangle we study, one can 
show that there are only three possible areas : A = l/v2, (v2 + l)/2 or (v2 — l)/2 (we take the 
hypotenuse of the triangle as unit length). We emphasize however that the formulae obtained in the 
present paper are quite general. 

The numerical spectrum was obtained by expanding the wave function near the angle vr/S in 
"partial waves" which are Bessel functions times a sinusoidal function of the angle : ip{r, 9) = 
^mmax Jg^(/j7-) gin(8m0). This automatically fulfills the Dirichlet conditions on the two faces of the 
billiard that meet at the corner with opening angle vr/S. The boundary condition on the remaining 
face is enforced in a manner identical to the improved point matching method presented in |22]. This 
results in a secular equation whose solutions are the eigenlevels of the system. We have tested the 
numerical stability of our procedure by varying the number rrimax of partial waves. We have computed 
the first 20 000 eigenlevels and we have checked that they were determined with an accuracy of the 
order of 1/1000 of the mean level spacing. 

The agreement with the numerically determined spectrum can be checked by studying the regu- 
larized Fourier transform on the level density : 



F{x) 



P 



7r{Aky 



dEp{E)e"'''e 



ikx -p{^^i^Y 



(8) 



In ^ krain and kmax are the lower and upper boundary of a window of the spectrum (typically kmm 
is the first eigen-level, and k^nax the 5000*^ one); kav = {kmax + kmin)/'^ and Afc = {kmax - kmin)/'^- 
/3 is a dimensionless regularizing parameter (typically /3 = 5). If kmin = 0, kmax — > +oo and /? = 0, 
F{x) is a series of delta peaks centered on the lengths of the classical and diffractive periodic orbits. 

The comparison of the result of Eq. (0) with the numerical datas is shown in Fig. ^ for the family 
of Fig. ^. The agreement is excellent. Note that in this figure (and in the followings of the same type) 
we compare different estimates for |F(x)|, but we also plot the modulus of the difference between 
the numerical F{x) and our analytical formula, which is a strong test of accuracy. Note also that 
for avoiding spurious sources of discrepancies with the numerical result we compute the integral (P) 
numerically even when we use an analytical expression for p{E) (this corresponds to what we still call 
the analytical F{x)). 



numerical 

S. CI. + corrections 




Figure 4: Comparison of the numerical 1-^(2;) | (dashed line) with the result from Eq. (0) (solid line). The two 
curves cannot be distinguished, this will also occur in all the following plots of the same type. We have taken 
here (3 = 5, kmin and kmax being respectively the first and the 5000"* level. The peak corresponds to the length 
of the family shown in Fig. ||, ^ = V2 (in all the text, the hypotenuse of the triangle is chosen as unit length). 
The modulus of the difference between the numerical and analytical values of F{x) is also plotted in this figure, 
but is barely seen on this scale (its largest value is 5 x lO""*^). The usual semi-classical contribution (with only 
the first term of the r.h.s. of (M)) gives instead an error of about 10~^. 



• We now concentrate on the second iterate of the family. Its contribution is less trivial than (^) 
and more generic ; hence we will present the computation in some details. One has here to consider 
double diffraction near the optical boundary. Eq. (^) is generalized to double (and in a similar fashion 
to multiple) diffraction : 



G2ir,r-",E) 



+00 



ds Gi(s , r ', E) n, ■ VgG^{r, s, E) 



+ Gi(fo, r ', E) Vregie, 9') Go(f, fo, S) , 



(9) 



where G2 is the total (semi-classical plus diffractive) Green function. 

When unfolding the trajectory (as done for instance in Fig. ^) near the diffractive boundary of the 
family, one is lead to consider contributions such as presented in Fig. ^. On that figure the position of 
a point r in vicinity of the diffractive trajectory on the boundary of the family is defined by coordinates 
X and y. x is a coordinate along the orbit {0 < x < i) and y a is transverse coordinate. 




Figure 5: Schematic representation of the different contributions to the Green function G2 near a diffractive 
periodic orbit on the boundary of the second iterate of a family. In this plot (and in the foUowings of the 
same type) the trajectory is represented after unfolding and the shaded areas are zones forbidden by classical 
mechanics. The path represented by a solid line going from r to si, S2 and r contributes to the leading KirchhofF 
term in the first integral of the r.h.s. of (||). The dashed path contributes at next order (it has one "Keller 
bounce" at tq with diffraction coefficient I?reg)- Its contribution to the Green function corresponds to the second 
term in (|lj). 

The leading order contribution of G2 to the level density is the usual contribution of the second 
repetition of a family of periodic orbit. It is obtained by simply making the approximation G2 ~ Go 
and it is of the form : 



p{E) 



A 



1 



2iT y/27rkL 



cos{kL — it/4:) . 



Here and in the following of this section L = n£ is the total length of the trajectory, 
primitive length, and n is the repetition number (here n = 2). 

The diffractive corrections to ( p!o| ) are included into p{E) through the following trace : 



(10) 



is the 



piE) 



TT 



Im (fr {G2{f,r,E)-Goif,f,E)) 



(11) 
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From the expression (^) of G2 one obtains the first order contribution to G2(f , f , E) — Go{f , f , E) 
under the form : 



^^ikL + 3i7r/4 



2{2Tr)y^y/xi{e-x) 



Jo 



.k f {si-yf , {s2-s,f , {y-S2)'' 
Q{y) I dsi / ds2e 2 I ^ ^ ^"^ 



—00 J —00 



• (12) 



This contribution corresponds to the main Kirchhoff term in (P) (i.e. to the sohd path from r to 
r in Fig. (^)) from which the semi-classical contribution has been removed when it exists, i.e. when 
y > (this semi-classical contribution has been written as the double sum J_ °° dsi f_ °° ds2 in ([l2|)). 



In the expression (|T2|) one has made the hypothesis |y — si|, |si — S2I, |y — S2I <C x,i. Thus, for instance, 
|s2 — si| ~ £ -|- (s2 — si)^/(2^). The above expression has to be inserted into Eq. (0), i.e. integrated 
transverse to the orbit (along y) and longitudinally (along x) . This is done in [Appendix B and the 
resulting contribution to the level density is : 



p{E) 



e 



SiT^k 



cos{kL) 



(13) 



This shows that the main diffractive correction to the contribution of the second iterate of a family 
is of order of the one of an isolated periodic orbit. Such non-generic contributions in vicinity of a 
family have already been studied in a slightly different context in Ref. |l( 



For a better agreement with numerical data, one needs to include also the next order correction 



to (13) in the level density. This is done by including mixed Kirchhoff-Keller contributions in the 
Green function (^), such as described by the path represented in Fig. (|5|) by a dashed line. Along that 
path, the first diffraction at fo is of Keller-type (involving a coefficient Vj-eg), and the second one of 
Kirchhoff type (with an integral along S2). One has two contributions, one for each possible location 
of Keller diffraction (a single one being shown in Fig. (|5|)). The relevant contribution to G2 are now 
of the form : 



T) J,kL + iir/A — if(77r/2 
AVk{2^f/^yJxl{l-x)~ 



-00 ,-fcJ (n-i/)^ I "1 I v'^ 



dsi e 



+ I ds2e I- 




(14) 



The integral of this expression is computed in Appendix B| (Eq. (|B10[) ). It yields the next correction 
to (O) which is of the form : 



p{E) 



V. 



Teg 



2Tik VStt/cL 



cos{kL — fd7r/2 — 37r/4) 



(15) 
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In ( p^ ) frf is the Maslov index of the orbit corresponding to the dashed path in Fig. 0. If o" is the 
relevant index near the optical boundary considered (i.e. the one for which 'Dfj,ri in Eq. (y) diverges), 
one can show that exp{— iz/(i7r/2} = a. 



There is a last correction to ( [l5|) , purely of Keller- type, giving a contribution of order 0{k~'^), but 
the contributions (10), (13) and ([l5|) already give a very good description of the Fourier transform of 
the spectrum. This can be checked in Fig. ^ for the second iterate of the family drawn in Fig. ^ 



numerical 

S. CI. + corrections 

S. CI. 



fa, 




Figure 6: Same as Fig. 4|fo r the second iterate of the family. The solid curve (labelled "S. CI. + corrections") 
corresponds to Eqs. (^ 13,1^). In this plot and in the foUowings of the same type, we represent with a thin 



dashed line (denoted "S. CI.") the usual semiclassical result without diffractive corrections (which corresponds 
here to Eq. (|o|) alone). The modulus of the difference between the numerical and analytical F{x) is also plotted 
in this figure, but is barely seen on this scale ; its largest value is 5 x 10"'*, whereas the usual semi-classical 
approach gives an error of 1.4 x 10^^. 



A simple remark is in order here : Eq. (|To|) is actually the first term of an expansion in k~^ (or 
equivalently in K). The magnitude of the next correction can be estimated by the following argument : 
the exact Green function in an infinite wedge can be expressed in terms of a Hankel function (with 
diffractive corrections unimportant for the present discussion) and this yields instead of Eq. (^) to 
something like 



p{E) 



^UkL) 



A 



1 



2-K V27rA;L 



cos{kL - 7r/4) + 



i A/{8L^) 
2TTk VSvr/cL 



sm.{kL - 'k/A) + 



(16) 



Equation (16) is the exact contribution of a family in an integrable polygon. Therefore there are 
non-diffractive corrections to the leading order ([l^) of the trace formula which are of same order as 
(Il5|). To V^eg in (15) one simply adds a factor 0{A/L?') (see the last term of the r.h.s. of ([TBD). 
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In all our numerical checks this correction appeared to be negligible. By comparing the diffractive 
corrections with (|l6| ) one can note that : (i) the first diffractive term ([l3|) is the leading \/% correction 
in the trace formula and (ii) for long orbits (or large repetition number) the term Ajl? in ( |T^ ) will 
be dominated by V^eq in ([l5|). 



• The determination of the contribution of the next iterates of a primitive family of periodic orbits 
with a diffractive boundary is patterned on the above derivation. We just state here the results. 

The main contribution is the generic semi-classical one, of the form (10). The first correction is of 
a type similar to the contribution to the trace formula of an isolated periodic orbit : 



p(E) < -CnCos{kL) , 

vrfc 



(17) 



where i is the primitive length, n is the repetition number, L = ni and Cn is a dimensionless parameter 
given by the formula C„ = (l/Svr) X]g=i [^(^ ~ q)]~^ ■ We show how to compute it in some special 
cases in [Appendix B| and in general in [Appendix D| . Its first values are Ci = 0, C2 = 1/(8 vr), 
C3 = l/(47rv2) ... and it has the limiting value Coo = 1/8- 

The next correction to (|l7) is of the form (|l5|). This is proven in special cases in Appendix B[ (Eqs. 
(BlCjBll)) and in general in Appendix C . 



We have tested the excellent agreement of contributions ([10| , [17[ , [15[) with the numerical spectrum. 
We illustrate this for the fifth iterate of the family shown in Fig. |^. This family is particular in the 
sense that one of its boundaries is formed by an isolated orbit which has an extra bounce compared 
to the family ; it lies along the lower edge of the triangles in Fig. |^. The contribution of such an 
isolated orbit has been known for some time |24,25| and is taken into account in the comparison with 
numerical results displayed in Fig. |^. The other boundary is a diffractive orbit of the type we are 
interested in. Its contribution to the level density is described by Eqs. (17,15). 




Figure 7: Representation of a family of periodic orbits in the triangle (7r/2, 7r/8,37r/8) by the technique of 
unfolding. The family has a length i = 2 cos 7r/8. Its area is shaded and the diffractive point on its boundary 
is marked with black dots. 



• The following diffractive corrections to the contributions (10,ni5) to the n iterate of a family 
correspond to orbits having n^ Kirchhoff diffractions and Ug Keller ones, with n^ + Ug = n. They 

rip + l . . 

yield corrections of order 0{k 2~) compared to the leading term ([10[) . One can show that their 
contribution to the level density is of the form : 
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■ numerical 
S. CI. + corrections 
S. CI. 




■ numerical 
S. CI. + corrections 
S. CI. 



9.240 9.250 

X 




Figure 8: Same as Fig. || for the fifth iterate of the family shown in Figure |^ (L ~ 9.239). The analytical 
result corresponds to Eqs. (E0^,n7l|l3) with n = 5. The modulus of the difference between the numerical and 
analytical F{x) is represented by the shaded area. It is important on the left plot, due orbits in vicinity of 
the peak that are not taken into account. This overlapping of peaks can be suppressed by increasing kmax (in 
order to decrease the width of the peaks in |F(a;)|). This is done on the right plot which is drawn for kmax 
corresponding to the 20 000'^ level. Then, the largest discrepancy with the numerical result is 4 x 10""* whereas 
the error when employing the usual semiclassical approach is 100 times larger. 



P{E) 



D 



reg 



27r/c \^8irki 



si""^ cos{kL - 3ng7r/4 - z/dvr/2) . 



(18) 



The Maslov index Ud in ([iq ) is related to the index a of the optical boundary considered (i.e. the 
one for which 'Da-,r) in Eq. (^) diverges) by exp{— zfrf7r/2} = o""". Bn ^ is a dimensionless coefficient. 






(1) 



n 



-1/2 



>(2) 



^n-1 - 



In agreement with (fL5|)), Bn = X]n=i 9 



3/2, 



n 



q) ^' ^ and the general form is : 



bic^^ = e n (^^^+1 - ^^ 



-3/2 



(19) 



with the convention Qug+i = n + qi. The sum is extended over all possible sets of Ug integers {qi}i<i<ng 
with 1 < qi < q2 < ■■■ < Qng < n. 

4 A diffractive orbit on the frontier of both a family and of the 
bilUard 

In the previous Section we have studied the case that a diffractive periodic orbit lies on an optical 
boundary corresponding to the frontier of a family. There is a special configuration where such a 
diffractive orbit lies on two optical boundaries. From Fig. || one sees that two optical boundaries meet 
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only on the edges defining the diffractive corner (when or ^' = or 7). Hence, in that case part of 
the diffractive trajectory crawls along the frontier of the billiard. This happens for instance for the 
diffractive trajectory on the boundary of the family shown on Fig. ^. 

Although the diffractive periodic orbit considered bounds the first iterate of a family, it is already 
doubly diffractive. For incorporating such a configuration into the trace formula, one can still use 
Eqs. (|5|) and @, but the semi-classical Green function to be incorporated in that formula has two 
contributions : one from the "direct" path (we call this path direct, but it may have bounces that 
are not materialized by the process of unfolding) and one from a path bouncing on the frontier of the 



billiard which is also the frontier of the family (see Fig. 10) 





Figure 9: Representation of a family of periodic 
orbits in the triangle {tt/2, tt/S, Stt/S) by the tech- 
nique of unfolding. The family has a length L = 
(4 + 2\/2)1/2 ~ 2.613. Its area is shaded and the 
diffractive points on its boundary are marked with 
a blaek spot. The boundary of the family party 
coincides with the frontier of the billiard. 



Figure 10: Schematic representation of the dif- 
ferent contributions to the Green function near a 
diffractive periodic orbit on the boundary of a fam- 
ily. The boundary of the family coincides with the 
one of the billiard along a segment of length £2- In 
this ease, the main Kirehhoff term in (O) contains 
two paths which are represented with solid lines of 
the figure. 



In that figure one has represented a configuration where a point in vicinity of the diffractive periodic 
orbit lies along the part of the boundary of the family which does not coincide with one frontier of 
the billiard (we denote by ii the length of this part, and by £2 the length of the part along a frontier 
of the billiard, ii + £2 = L). Then, the main Kirehhoff contribution to G2 is of the form : 



^/Z QikL+Sin/i 



2(27r)3/V^^2(^i-x) 



00 ^+00 ^ 

dsi / ds2 e 
Jo 



+ 00 ,-fe J ("i-yf (^2+''l)^ I fa-''2)^ 



dsi / ds2 e 
Jo 



(20) 



The second contribution in ( [20| ) is obtained from the first one by the method of images. It 
corresponds in Fig. |l^ to the path going from si to 5*2 with one bounce on the frontier of the billiard. 
If the point r of Fig. |lO| lies along the part of the orbit coinciding with the frontier of the billiard. 
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then the main Kirchhoff contribution to G2 is a sum of four terms (this is detailed in [Appendix B , 
cf. Fig. 24). We will not give the explicit computation here (see Appendix B| ), but after transverse 
integration along y the final result for the first diffractive correction to the contribution of a family 
such as the one presented in Fig. |9| is : 



p{E) 



An'^k 



\Jhh + ^ arctg 




cos{kL) 



(21) 



The next diffractive corrections to (^) are of the order of a doubly diffractive Keller correction, 
and we do not include them in our description of the family. As seen in Fig. 11, contributions (Il( 



and (|2l|) already give an excellent description of the Fourier transform of the spectrum in vicinity of 
the length of the family drawn on Fig. (0). 







Figure 11: Comparison between the numerical evaluation of F{x) with the result of Eqs. ( |l0[|2l| ) for x close to 
the length of the family shown in Fig. g. The modulus of the difference is also plotted on the figure, but cannot 
bee seen (it is lower than 10""^). 



5 A diffractive orbit jumping from the boundary of a family to the 
boundary of an other one 

In the billiard we consider, an interesting combination of orbits occurs. It is formed by the gathering 
of two diffractive orbits, each being on the boundary of a family, and where the total diffractive orbit 
is on the optical boundary, although the two families have no overlap. An example of such a case is 



given in Fig. 12 



As seen in Figure |12|, although the diffractive orbit lies on the optical boundary, there is no allowed 
classical trajectory nearby. This type of orbit might be a particularity of Veech billiards, but it is 
nevertheless interesting to describe its contribution to the level density. The schematic representation 



of the neighborhood of the orbit is displayed in Fig. 13. 
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Figure 12: Representation of a difFractive orbit 
that jumps from the boundary of a family to the 
boundary of an other one. The area occupied by 
each family is shaded. The first family transports 
triangle A to position _B, and the second one trans- 
ports triangle B to position C. The difFractive 
orbit is the thick solid line with its three diffrac- 
tive points marked by black dots. It has a length 
Ld = (10 + 7^2)1/2-4.461. 



Figure 13: Schematic representation of the con- 
tributions to the Green function in vicinity of a 
difFractive orbit that jumps from the boundary of a 
family to the boundary of another one (such as the 
one shown in Fig. |l^). ti corresponds to the length 
of the family that maps triangle B onto triangle C 
in Fig. [l^ and (.2 + ^3 to the length of the family 
that maps triangle A onto B. For this last family, 
(.2 is the length of the part of its boundary that lies 
along the frontier of the billiard. 



In the cases we have studied, the problem is complicated by the fact that one of the families 
considered has a boundary that partially coincides with the frontier of the billiard, i.e. is of the type 
studied in the previous section. Hence, there are three relevant lengths along the orbit we consider : 
ii is the length of one of the families, £2 + ^3 is the length of the second one, (.2 being the length that 
corresponds to the part of the boundary of the second family lying along the frontier of the billiard. 

The difFractive Green function to be considered here is of the the type previously studied in Section 
0, with an additional diffractive bounce. Hence, one defines a Green function G3 connected to G2 in 
the same manner as G2 is connected to Gi in Eq. (^. Due to the possible bounce along the frontier 
of the billiard that coincides with the boundary of the family, the semi-classical Green function to be 



incorporated in this formula has several contributions. This is illustrated in Fig. 13 where there are 
two possible paths for going from si to 5*2 (such a contribution was already present in Fig. 10). If 
the initial point r were lying near the frontier of the billiard (next to the part of the family of length 
£2), one would have four different paths : two for going from r to 5*2 and two for going from si to 
r . After transverse integration of these four contributions (along the variable y), one can verify that 
they lead to the same contribution as the ones displayed in Fig. fsl hence we will only present here 



the computation in the simpler case shown in Fig. 13 



For the configuration of Fig. 13, the Kirchhoff term in the expression of G3 reads : 
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{2ikf D^{x,h,h,^ 



oo f+oo rO 

dsi / ds2 / (isfie 



1 / af,2 / assi 
Jo J-oo 

+00 ^+00 /•O , 

ds^ / ds2 / ds.-? e 



(ai-y) I (g2~°l)" I ("3~"2) I (v-i's) 
X ^ I2 ~'~ ^3 "■" fi-i 



1 / as2 / as3' 
Jo J-00 



X "■" fo ^ fQ "■" fi -I 



(22) 



where the notation Dn is defined in [Appendix B| . In (P2|), the second term of the r.h.s. is obtained 
from the first one by the method of image. It corresponds to a path going from si to S2 with a specular 



bounce off the frontier of the bihiard (cf. Fig. 13). Note that there is no classical path contributing 
to (p^) : it is clear from Fig. 13 that there exits no classical trajectory from f to r . The transverse 



and longitudinal integrations are done in a manner similar to the one exposed in Appendix B for the 
similar case of an orbit which boundary coincides with the frontier of the billiard (cf. Section ^. The 
contribution of the diffractive orbit to the level density is : 



p{E) 



L„ 



Sn^k 



cos{kLii) < arctg 



h+h 



arctg 



is + h 



+ arctg 



h+i2 



VWT+h) , vW^+h) Vh{h + i3) 



where L 



d 



Ld Ld 

-1 + ^2 + ^3 is the length of the diffractive orbit. 



L„ 



(23) 



In order to have a good description of the contribution of the orbit we are considering here, one 
needs (as in Section 0) to incorporate next order corrections, i.e. mixed Kirchhoff-Keller terms. This 
corresponds in Fig. |l^ to the path with one Keller bounce on the apex which is not on the frontier of 
the billiard (dashed line) (Keller bounces on the other apexes contribute to higher order). We do not 
detail the computation here and just present the resulting correction to (P3|) . It is of the form : 



P{E) 



Ld V, 



reg 



2TT^k ^SirkLd 




/ t3"l 

arctg W cos{kLd — v^t^ 12 — 37r/4) . 



2-^d 



(24) 



As one can see in Fig. 14, Eqs. (p3|,24) give an excellent account of the contribution to the level 
density of the orbit shown on Fig. |12|. 



6 A diffractive orbit bouncing between the upper and lower boun- 
dary of a family 

Up to now we have only considered diffractive orbits which were lying exactly on the optical boun- 
dary. Other types of diffractive orbit occur which do not stand right on the optical boundary, but 
close enough to prevent their description by the geometrical theory of diffraction. Such an orbit is 
represented in Fig. |^. 
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Figure 14: Same as Fig. for tlie orbits shown of Fig. n3. The sohd line corresponds to the results of 
Eqs. (p3,E3). The shaded area at the bottom corresponds to the modulus of the difference between the numerical 
and the analytical result (which is inferior to 2 x lO""*). Note that the standard semi-classical approach 
completely misses this peak in |i^(x)|. 

On this Figure, the upper left triangle and the lower right one are connected by a family. For 
legibility we do not represent it and its area. We only represent the diffractive orbit on its boundary 
(the straight line between two black dots). This orbit is singly diffractive and its contribution corrects 
the one of the family as in Eq. (^). There is a diffractive orbit nearby, not exactly on the optical 
boundary, but very close to being part of the family : it starts and ends at the same point as the 
diffractive orbit on the boundary of the family, but it has an extra diffractive bounce in between (see 
Figure [l5| ). This is the type of orbits we aim to describe in this section. Its diffraction coefficient 
does not exactly diverge, but one of the contributions ^o-.r; of (|3|) is large and does not allow a proper 
description of the diffractive Green function by means of Eq. (g). For simplicity we will denote this 
part as the "divergent part" (the remaining being the "regular part"). 

The configuration we just described is of the type represented schematically on Fig. ^. The 
projection of the upper diffractive corner onto the family separates it in two parts of lengths ii and 
^2 (^1 + ^2 = L). Typically, the upper wedge represented in that figure is the upper boundary of the 
family. In that case, if (" denotes the distance from the upper wedge to the lower extend of the family, 
the area occupied by the family is simply A = C L {L being the length of the family, see Figure ^). 

In this configuration the leading term in the Green function is obtained from the explicit expression 
of G2 and reads : 



G2(f,f,^) 



^^ikL + 3i7r/A 



2(27r)3/2^a;^2(^i-a;) 



dsi 



+00 
dso e 



;fej in-v)'^ I (a2-ai)^ , (V-S2)^ 



(25) 



The above expression integrated transversely (along y) and longitudinally (along x) gives the 
contribution of the family and of its corrections to the level density. The result reads (the relevant 



integrals are given in Appendix B) 
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Figure 15: A difFractive orbit on tiie boundary of a 
family (straight line between two black dots). The 
family is not represented, but it connects the upper 
left triangle to the lower right one. The other orbit 
shown is typical of those studied in this section. It 
is a doubly difFractive orbit close to the family, it 
is represented by the segmented line between three 
black dots. 



Figure 16: Graphical representation of different 
diffractive contributions for a family limited by two 
opposite wedges. The leading term in m) cor- 
responds to the solid line. The dashed path is 
one of the next order corrections, involving one 
Keller bounce (with a coefficient T>reg) on one of the 
apexes. C is the distance between the two apexes 
measured transverse to the direction of the family. 



-Im I d^rG2{f,f,E) 



+ 






A 



2tt V27rA;L 



cos{kL — it/4:) 



cos{kLd) - 2 V^^kA Re le^'^^t-i^/^K 



(26) 



In (|26D, K is the modified Fresnel function defined in [Appendix A| . We have denoted by L^ the 



length of the doubly diffractive orbit going from the upper corner to the lower one (L^ = v^i + C^ + 
Y^2 + C^)' ^y ^ ^^^ length difference L^ — L, and have made the approximation A ~ ("^L/(2^i^2)- 
The first term of the r.h.s. of (^) is the usual contribution of a family. The second is a diffractive 
correction. 



Two remarks are in order here 



• It may happen that the upper corner of Fig. |T^ does not provide the upper boundary of the family 
because the family meets an other non-diffractive boundary between the two diffractive corners. This 



is the case presented in Fig. 15 ; the family does not occupy all the width between the two diffractive 
corners : it meets first a non-diffractive 7r/2 corner. In this case formula (|2q ) remains valid, but 

A = CL/2. 



• Secondly, it is interesting to check what is the behavior of Eq. (26) when the two wedges are far 
apart, i.e. in the limit \^kA S> 1. By using the asymptotic expansion (|A3| ) of the modified Fresnel 
function one obtains : 
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(27) 



In (27), to the usual contribution of a family is added a term which can be matched with a 
contribution such as (Q) with two diffractive bounces, provided some approximations are made. The 
term 2iii2/{L(iA) stands were one would expect a product of two coefficients P. Indeed, one can show 
that this term corresponds to the product of the two divergent parts ^o-.r; near the optical boundary. 
But it is not of the form (^) which is the only acceptable in the limit where Eq. ( p7[ ) has been written 
(i.e. far from the optical boundary). This is a well known drawback of Kirchhoff's approximation 
already discussed in Section |2.2| . It can be cured relatively easily : if the optical boundary close to the 
diffractive orbit is characterized by the indices a and t], one has to multiply the second terms of the 
r.h.s. of dH) by the factor (|a^,^| V^^^f/2 and to express A as ja^^^p^i^a/^d (instead of Ld-L). The 
term ao-,r] appearing in these expressions is defined in Eq. (1A4D . The upper index d is meant to remind 
that Oo-,r; and Po-.r? have to bee evaluated on the diffractive periodic orbit of length Ld- This procedure 
allows recovery of the correct limit in (p^). Moreover it does not affect (26) when the diffractive 



orbit is close to the family (i.e. in the limit vfcA <^ 1) since in this limit (|a^„| 2?^„) /2 ~ 1 and 
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Figure 17: Same as Fig. Hfor the orbits shown on Fig. UB. The solid line corresponds to Eqs. (Eq,E8[). Here one 
has L = (1 + V2)V6 ^ 5.9136 and Ld = (10 + 3^2)1/2 + (6 - V2)^/^ ~ 5.9154. The shaded area hardly seen at 
the bottom of the plot corresponds to the modulus of the difference between the numerical and the analytical 
result (which is inferior to 1.4 x 10"'^). The pure semi-classical estimate (Eq. dlQ)) gives an error of 4.6 x 10~^. 
Taking only into account the diffractive periodic orbit standing exactly on the boundary of the family (as in 
Eq. (^) gives an even larger error (5.3 x 10^^). 

Eq. ( p6[ ) is not the final contribution from the configuration represented in Figure 16. This is clear 
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from (|2^) : far from the optical boundary the asymptotic evaluation of ( |26| ) only allows to recover 
the divergent part of the diffraction coefficient. Hence one has to include other terms, of mixed type 
Keller-Kirchhoff, as already encountered in Sections ^ and |5|. These involve the regular part Dreg 
of the diffraction coefficients. We must be careful though, that one has now two different diffraction 
coefficients : one for the orbit along the boundary of the family (we denote it by Pfeg) and an other one 
for the orbit bouncing from the lower wedge to the upper one (we denote it by D'^ ). The remaining 
contribution to p{E) of the configuration considered in this section is (we do not detail the derivation) : 

P^^^ ^ :^-^=cos{kL-u,n/2-37T/4) 
2TTk VSvrfcL 



T V^ 

vrA; ^/8^TkL 



{e^' 



2^kSTTk./hh "• ^ ' > ^ ' 

The term involving a modified Fresnel function in the above expression can be made uniform by 
a procedure similar to the one devised for Eq. (p6|). Note also that we have added in ( |2^ ) a doubly 
diffractive term of purely Keller type (last term of the r.h.s.). It is a small correction and such terms 
were neglected in the previous sections. We kept it here for consistency because far from the optical 
boundary, it is of same order as the second term of the r.h.s. of (|26|). 

The agreement with the numerical spectrum is here also excellent, as shown by Figure |l^ Note 
that the geometrical theory of diffraction - although yielding a non-divergent result - is completely 
inadequate in this case. It amounts here to treating the isolated diffractive orbit as truly isolated from 
the family : hence to describing the family of periodic orbits in the usual way (i.e. using Eq. (p!o|)) and 
including a correction of type (Q) describing the contribution of the doubly diffractive orbit bouncing 
between the upper and lower boundary of the family. This procedure gives an error of 9.4 x 10~^ in 



Fig. 17. 



7 A diffractive orbit near an isolated one 

In this Section we will study, as in the previous one, a diffractive orbit standing not exactly on the 
optical boundary, but close to an allowed 'periodic orbit. Here we consider the case that the nearby 
orbit is an isolated one (and not part of a family as in the previous section) . Such a configuration has 
already been studied in Ref. [^, and we will here re-derive the result in a simpler manner (but with 
less generality). 



A typical occurrence of the situation we are interested in is shown on Fig. 18. The isolated orbit 
we consider is the third iterate of the shortest classical periodic orbit of the system. It has a length 
L = 3/^2 ~ 2.121. The nearby singly diffractive orbit has a length L^ = (6 - \/2)^/2 ~ 2.141. 



The different contributions to the Green function Gi are illustrated in Fig. 19. Note that for the 



phase-space coordinate transverse to the direction of an orbit, a bounce on a straight segment leads to 
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Figure 18: A diffractive orbit (straigiit line be- 
tween two black points) near an isolated orbit 
(straight line connecting two corners with opening 
angle 7r/2). 



Figure 19: In this figure the isolated periodic orbit 
connects the two open circles. A nearby path from 
r to r is represented. There is an inversion along 
the periodic orbit for the transverse coordinate y 
(see the text). The diffractive periodic orbit goes 
from one open circle to the apex at tq and to the 
other open circle. It has a length Ld — L + 2('^/L. 



an inversion. Hence, in a polygonal enclosure, the transverse mapping near a periodic orbit is either 
an inversion (for an odd number of bounces) and the orbit is then isolated, or the identity (for an 
even number of bounces) and the orbit is then part of a family. This is the reason for the inversion on 



Fig. 19 of point r with respect to the axis of the isolated periodic orbit after the process of unfolding. 
In the Figure, C is the distance from the diffractive apex to the periodic orbit. 

In the case of interest here, the Kirchhoff part of the total Green function Gi is (from Eq. (P) and 
Fig. [ID : 



JkL — ii'TT/2 



47ry^x(L 



dsi exp <j i— 



X L — X 



(29) 



where z^ is the Maslov index of the isolated orbit (exp{ii^7r/2} = — 1). If i^d is the one of the diffractive 
orbit and if a characterizes the nearby optical boundary, one has exp{ii'[i7r/2} = aexp{iuTT/2}. Once 
Go has been removed, the above expression yields - after transverse and longitudinal integration - the 
main contribution of the diffractive orbit to the level density. There is also a corrective term containing 
the regular part of the diffraction coefficient. Altogether one obtains the following contribution : 
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2irk^/8TTkL2 



cos{kLd - Vd-n/2 - 37r/4) , (30) 



where A = L^i — L c^:^ 2C,'^ /L. As in the previous section we have used a representation of the Green 
function based on Kirchhoff 's approximation which does not yield a uniform formula : Eq. ( pO[ ) does 
not permit recovery of the result of the geometrical theory of diffraction far from the optical boundary. 
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i.e. when the isolated and the diffractive orbit are far apart. As in Section |^, one can easily remedy 
this deficiency. If the optical boundary to which the diffractive orbit is close is characterized by the 
indices a and ??, one multiplies the first term of the r.h.s. of pO| ) by o"(Lrf/L)|ao-^^| V^j^ri/V^ and replaces 
vA in the argument of the modified Fresnel function by |ao-^,j| \/L2/2. This procedure does not affect 
( |30[ ) in the limit that the diffractive and isolated periodic orbits are close and it allows recovery of the 
result of the geometrical theory of diffraction when these two orbits are well separated. 






numerical 

S. CI. + corrections 




Figure 20: Same as Fig. ^ for the orbits shown on Fig . [l8| . The solid line corresponds to the contribution of the 
isolated diffractive orbit (which has a length L = 3/V2c± 2.121) plus the contribution of the nearby diffractive 
orbit {Ld = (6 — v^)"'^' ^ — 2.141). The shaded area at the bottom of the plot corresponds to the modulus of the 
difference between the numerical and the analytical result (which is inferior to 5 x 10^^). We have used here 
the uniformized version of Eq. (30) (see the text), the use of the plain formula gives twice a larger discrepancy 
with numerical datas. 



The comparison of formula (BO) with the numerical result is very good, as shown in Fig. 20. Note 



that the geometrical theory of diffraction is not totally inadequate here (as it was in the previous 
section). It gives an error only four times larger than our approach. The reason is that the classical 
and diffractive orbits considered here are not very close to each other. Of course, the distance between 
two orbits must be measured relatively to the wave length. As a result the accuracy of the geometrical 
theory of diffraction depends on the window of the spectrum chosen for evaluating F{x). For instance, 
evaluating F{x) keeping only the first 500 levels (instead of the first 5000 levels as in Fig. |20|) gives 
for the geometrical theory of diffraction an error 10 times larger than our approach. 



8 A rectangular billiard with a flux line 

In this Section we depart from the previous examples which treat corner diffraction, and we consider 
instead diffraction by a fiux line. We consider a rectangular billiard (with sides of length a and b) 
with a flux line located at point fg inside the billiard (cf. Fig. pi|a). 

We will not restart here a detailed study of a large number of different cases of diffraction in the 
system (as was done in Sees, gtol^for a triangular billiard). First, because Aharonov-Bohm diffraction 
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is in a sense simpler than corner diffraction and leads to fewer exceptional cases ; second because we 
chose this example merely to illustrate the flexibility of Kirchhoff 's approach devises in Sec. ||. We 
will show that Eq. (^) permits to tackle the problem of multiple forward Aharonov-Bohm scattering. 



(a) 



i 


i 




' 



( b 




Figure 21: Part (a) : the rectangular billiard and the two bouncing ball orbits (the vertical and the horizontal 
one). The flux line is located at point tq niarked by a white circle. Part (b) ; representation of the typical 
contribution to the Green function (g) in vicinity of the doubly diffractive orbit belonging to one of the bouncing 
ball families. The path going from r, to si, 5*2 and back to r represented in the figure accumulates a phase 
exp(— 2i7rQ;)) (see the text). 

This problem is encountered for instance when evaluating the contribution to the trace formula of 
the two families drawn in Fig. pi|a. For each of these families the periodic orbit that encounters the 
point vq twice on its way gives a doubly diffractive contribution. The schematic contribution to (0) for 



a nearby closed path is illustrate in Fig. 21b. In this Figure there is a reflection on the frontier of the 
billiard between the two flux lines and this has the effect of changing the sign of a on the second flux 
line (equivalently one could keep the same a and change the orientation of the axis (ro,,S2))- From 
Eq. (^) the diffractive Green function of the problem is written : 



G2dif,f,E) 



2TT^8TTXe2{h -X) 



00 r+00 



dsidso e 



00 J —00 



gi7ra{Sgn(si)-Sgn(s2)} _ -^ 



(31) 



where i is the length of the periodic orbit. The flux line (encountered twice) separates the orbit in 
three parts having lengths denoted by x, £2 and ii — x in ( |3ll) and Fig. ^l| b {ii + £2 = (-)■ Transverse 
integration yields : 



+00 
dy G2d{r,r,E) 

-00 



Vhh iM + i7r/2 



2-nU 



{cos(27rQ) — 1} 



(32) 
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and this gives a contribution to the level density : 



p{E) 



VW2 



vr 



sin (vra) cos{M) . 



(33) 



We check in Fig. 22 the very good agreement with the Fourier transform of the spectrum in vicinity 
of the length of the families drawn on Fig. |2l|a. In this Figure, the numerical F{x) is computed using 
Eq. (P) with /3 = 5, kmin and kmax being respectively the first and the 1400 level. In the numerical 
computation we took a = 1/2 because in this case the diffractive effects on the level density are at 
maximum (see Eq. (|33|)). The shaded area hardly seen at the bottom of the plot is the modulus of 
the difference between the numerical and analytical F{x). For obtaining the excellent agreement of 
Fig. ^ we have taken into account classical isolated boundary orbits (of the type already encountered 
for the family drawn in Fig. ^) and simple nearby diffractive periodic orbits which can be treated 
within the geometrical theory of diffraction (the relevant formulae are given in [^ ) . Not taking into 
account the diffractive contribution (|3^) would give a much larger error which is represented by a thin 
dashed line. 
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Figure 22: Comparison of the numerical |i^(x)| (dashed line) with the result from Eqs. dlfl,p3) (solid line). The 
large peaks correspond to the lengths of the families shown in Fig. Gm. In our computations we have taken the 
sides of the rectangle to be a = 4 and b — it and the bouncing ball families have then a length i = % and 2-k. 
The shaded line barely seen at the bottom is the modulus of the difference between the numerical and analytical 
F{x). We also show as a thin dashed line the value of this difference when not including the contribution (p3). 



9 Conclusion 

In this paper we have studied diffractive corrections to the semiclassical trace formula for the level 
density of polygonal billiards. Special care has been devoted to the treatment of diffractive periodic 
orbits lying on (or in vicinity of) the optical boundary, i.e. on the verge of being allowed by classical 
mechanics. In particular we derived a systematic expansion for the corner-diffractive corrections to 
the n^^ iterate of a family of periodic orbits. 
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The method employed (based on approximation (^) allows to treat a rich variety of different cases 
with great precision (Sees. |2|to|3)- This method is easily extended to similar diffraction problems. In 
particular, our approach to the diffractive correction of the n iterate of a family allows to treat the 
non-trivial problem of multiple forward Aharonov-Bohm diffusion (Sec. P). 

The main purpose of our study was to establish the basis of a trace formula in pseudo-integrable 
systems, with contributions from diffractive orbits. It seems that these diffractive corrections are 
responsible for particular forms of spectral statistics observed in many such models [^I|. Further 
investigations will elucidate this relationship. 
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Appendix A 

In this Appendix we derive Kirchhoff-like formulae for the Green function in the cases of corner and 
flux line diffraction (Eqs. (^,^)). Here we compare Eqs. (|^) and (^ with the exact diffraction in the 



free plane : diffraction by an infinite wedge (in Appendix A.l ) and by a flux line in the plane (in 



[Appendix A. 2]) . Adding boundaries to the problem (e.g. putting the flux line in a billiard) amounts 
through the method of images - to add other sources of diffraction. In this case we describe multiple 
diffraction by using a natural generalization of Eqs. (^J^ (see e.g. Eq. (^)). 

Appendix A.l Corner diffraction 

A uniform approximation for diffraction on a single corner has been first given by Pauli |21]. The 



subject has been studied in detail in the late 1960s and in the 1970s. We state here the result of one 
possible method of "uniformization" (detailed derivation and references can be found in ||lT 



Gid{r, f', E) = - — ^==^ 2^ \aa,n\ 1^a,r, K 
4 Vvr/cLd ^^^^ 



^o",»?l 



kW 



(Al) 



where £' and i are the lengths of the classical trajectories from f to rg and from tq to r (L^ = i' + i). 
X is a modified Fresnel function which is defined by 



K{z) = -= exp {-iz2 _ ivr/4} / dte'' = —— erfc (e"^"/^ z) , (A2) 

and which has the following limiting properties : K{0) = 1/2 and 

g«7r/4 i 3 TT Stt 

-^(^) ~ ^^/^ (1- ^-^•••) when \z\^+co with - - < arg (z) < — . (A3) 
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In (|A1|), Co-^jj is a kind of a measure of the angular distance from the trajectory to the optical 
boundary. On the optical boundary characterized by a and rj, a^^r} = 0. Far from the optical boundary 
its precise value is irrelevant since one can use Keller's approximation (which corresponds to keeping 
only the first term in expansion ([A3|) ) . In the transition region one has to use a specific form of aa^-q, 
which characterizes the type of uniform approximation chosen. We take here (see |17|) 



Oo-,r; = V2 COs( 



"-(T.r; 7) 



where 



rifj^ri = iimt 
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gz, 



(A4) 



nint denoting the nearest integer and (pa = 9' — (t9, where 9' (9) is the incoming (outgoing) angle of 
the diffr active trajectory with the boundary. 

In the following of this Appendix, we will use the uniform approximation ( |A1| ) for justifying the 
approximation (^) which is valid in vicinity of an optical boundary. Let's consider that for one of 
the four couples of values of (a, rj) one is near the optical boundary. In the contribution of the three 
other terms in (|A1| ), the modified Fresnel function can be evaluated by keeping only the first term in 
expansion ( |A^ ) and this gives the second contribution in the r.h.s. of (|5[). 




Figure 23: Representation of a diffractive trajectory going from f' to r. Part (a) displays the notations of 
Eq. (A5) and part (b) of Eq. (A7). In the figure, riS(l)cr > (or equivalently y' /d' + y/d > 0) and the classical 
trajectory from r to r is allowed. 



one needs to make explicit computations, 
for an orbit near the optical boundary (the 



For obtaining the first contribution of the r.h.s. of ( 

The configuration we study has been represented on Fig. | 

optical boundary for trajectories issued from f' is the dashed line of Fig. [2^). Note (from Ref. p7| ) 
that the classical orbit on the optical boundary has properties depending on a and rj. If o" = 1 (—1) 
it has an even (odd) number of bounces near the corner ; if r/ = 1 (—1) it bounces first on the line 
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= 7 (0)- If one writes (pa- = (paO + ^(pa- (where (paO = '^na^nl — ijir is the value of (pa on the optical 
boundary), by examining the four different configurations, one can convince oneself geometrically on 
Fig. |2^ that the oriented angle between (rg — f ') and (r — vq) is rj 6(pa- If rj 6(pa < 0, there is no classical 
orbit from f' to r. If rj6(pa > 0, the classical orbit is allowed and it has a Maslov index ly such that 
exp{iz^7r/2} = a. Since one is near the optical boundary, the angle 5(pa is small and in (Al) one can 



make the approximations a^^n ^ ??<5</'o-/v2, ^o-,r; ^ 2(T77/(5</>o- (compare with the exact formulae ( |A4| ) 
and dsh). One has also £ + £' — \aa,-q\'^U' /Ld ~ L, where L is the length of the classical path from r' 



to r. Altogether one obtains from (lAll) 



GUr,r',E) ^ sgn(r?5</),) / e^* dt 

2Tr\/2kL J^k{La-L) 

+ Go{r,fo,E)Vreg{e,e')Go{ro,r',E), (A5) 



and 



Gi(f, r ', ii;) = Girf(f, r ', £;) + e(r? 5</.,) Go{r, r ', S) , (A6) 

where is the Heaviside function. 

We will now show that this expression matches Eq. (|5|). For that purpose we will explicitly 
evaluate Eq. (||) by choosing a particular axis of coordinates, shown in Fig. p3|(b). In that figure we 
have chosen the locus of points s such that the distance from f' and f to its perpendicular is small 
(this is consistent with the fact that the trajectory is near the optical boundary). If one denotes by 
d {d') the distance between tq and the projection of r (f') on the perpendicular and by y {y') the 



algebraic distance from r (f* ) to the perpendicular (see Fig. |2^(b)), one has : 



ik{d+d')-iu-K/2 r+oo 



2/ dsGo{s,f',E)ns-VgGo{r,s,E) ^ / ds - 

Jo 47r Jq 



2 



pikL—iviT 12 f+oo . n 



2-K\/2kL J-sgn{riS(f>^)y/k{La-L) 



e" dt , (A7) 



where one has made the change of variable t = \/k{d + d')/{2dd')[s — {dy' + yd')/{d + d')] and one 
has used the facts that Ld — L k, {dd' /2){y' /d' + y/d)'^/{d + d') and sgn {y'/d' + y/d) = sgn {rjScpa). 
This expression inserted into (|5|) is equivalent to expression ( [A6| ) for the Green function. Hence, 
starting from the uniform approximation (Al), we have proven the validity of Eq. (||) near the optical 
boundary. 

Appendix A. 2 Flux line diffraction 

In the case of diffraction by a flux line, a uniform solution has been worked out by Aharonov and 
Bohm [^]. The corresponding expression for the Green function, as given by Sieber ^^, reads : 
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ikLa+i{4>-<l>')/2+iTT/4 

Gi{f, r ', E) = Go(r, r', E) e^°('^-<^ ) + sin(a7r) ^^^^ K 



l2kU' 



cos 



, (A8) 



where (p and (/>' are here the angular coordinates of a polar system with origin on the flux line rg- They 
have to be chosen such that |(/>' — (p\ < vr. The other notations are identical to those of Eq. (Al). In 
the conflguration illustrated in Fig. ^ this means that if 6(p is the angle between r — tq and vq — f' {6(p 

(j)' = Sep — TT sgn {5(j)) . Note also that as in Ref . p6| we 



is the analogous of rjScpa of Fig. 23a) one has 4> — (j)' = 5 

restrict ourself to the case of non-singular behavior near the flux line (i.e. vanishing wave- functions). 



We will now follow the same procedure as in Appendix A.l and show that Eq. (|6D and ( [A8D are 
equivalent in the limit of small 54>. In this limit the geometrical theory of diffraction fails, and indeed 
it is well known that the Aharonov-Bohm scattering amplitude diverges in the forward direction. 
The computations are very similar and for simplicity we chose here the arbitrary locus of points s 
perpendicular to the optical boundary which is the line going from f ' to rg , then y' = (the notations 



are deflned in Fig. 22 b). In the semiclassical limit, Eq. (y) reads : 



Gi{f,f',E) 



Jk{d + d!) + iai 



+00 ik_ 



(s-yr 



(A9) 



Using the fact that L k, d + d! + {l/2)y^ /{d + d') and Ld ^ d + d' + y'^ /{2d) this can be rewritten 



as : 



Gi{ry,E) 



J{kL + a6(/)) 



it I AocK 



dte'' +e 



/fcA 

dte' 



TTy/8k{d + d') \ J-eVkA 

where A = L^ — L and e = sgn {6(j)) = sgn (y). Simple manipulations show that 



(AlO) 



Gi{r,f',E)^Go{f,f',E)e 



ia{S(f>-€n) , Sin(a7r) i(kLd+TT /4-eTr/2) 



+ 



\/27rkL(i 



K VkA 



and this expression matches (A8) when \6(f)\ « vr. 



(All) 



Appendix B 

In this Appendix we give some useful formulae that correspond to transverse integration of the different 
types of Green functions appearing in the main text. 

Let's first define F{x, y, si, . . . Sn) by 



F{x,y,si,... ,Sn) =exp 



ik 



{Sl - yf ^ {S2-Slf ^ ^ {Sn - Sn-lf ^ {v 



'il-X 



(Bl) 
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and Dn{x,ii,...in) by 



Dn{x 



JkL-3(n+l)in/4 



, ti, . . . , tn) 



i7rk)'^y/x{ii-x)i2...ln ' 



(B2) 



where L = ii + . . . + in and < x < ^i 



• In Section ^, for treating the first order diffractive correction to the n repetition of a family, 
one needs to compute the following integral (cf. e.g. Eq. (jT^) which after transverse integration yields 
the expression of l2{i,i)) : 



/+00 p+oo ^+00 

dy / ... / dsi...dsn F{x,y,si,...Sn) 
-00 Uo Jo 

/+00 /•+00 
... j dsi...dsn F{x,y,si,...Sn) .(B3) 
-00 J —00 

It is easy to see that the x dependence disappears in expression (p^ after integration over y. This 
is the reason why we did not include x in the list of arguments of /„. 



We obtain 



pikL+iTT/2 



AnkL 






^ikL+iTT/2 



SnkL 

pjfcL+J7r/2 

levr/cL 



{Vhih + is) + VU^3 + h) + VUh + h)} 



(B4) 
(B5) 



{Vh{i2 + h+h) + VUh+h+h) + ^^3(4 +^1 +^2) 



+ V^4h +i2+ is) + Vih + ^2)(4 + h) + Vih + ^4)(^2 + h) 



+ -Vil{i2+i3+i4)aT^Ctg 

vr 



i2i. 



2<^4 



2+^3 + 4) 



+ -V^2(^3 + -^4 + ^i)arctg 
vr 



isii 



Uiis + i^ + ii] 



+ -\/^3(4 + 4+^2)arctg 



t4<^2 
h{i4+il+i2) 



+ -y/hihTi^Tis) aict^ 



iiis 



^2(^1+^2 + ^3 



(B6) 
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and 



2/5(^1,^2,4,^4,4) = /4(4+4,4,4,4)+/4(4 + 4,4,4,4)+/4(4,4+4,4,4) 



+ /4(4,4,4+4,4)+/4(4,4,4,4 + 4)-/3(4+4 + 4,4,4) 



h{h,i2 + 4 + 4,4) - /3(4,4 + 4,4 + 4) - ^3(4, 4, 4 + 4 + 4) 



+ /2(4,4 + 4 + 4 + 4) 



(B7) 



Although this is not apparent in the above expression, exphcit computation shows that formula 
(p7|) is - as Eqs. ( p^ , p5| , p6|) - invariant under cyclic permutation of the indices. 

Expressions ( |B4| ) to ( [B7D greatly simplify when all the ^'s are equal. One obtains : 



I^^i^,,,i) = ^eikL + i7r/2 

rv 



(B8) 



with Ci = 0, C2 = 1/{8tt), Cs = 1/(W2), C4 = (1 + 4/^)/(16 7r), Cg = (1 + ^/3)/(8 7r). A 
general formula for C„ is given in Appendix D| . From ( |B8| ) the contribution (O) to the level density 
follows directly. 



The next order correction to ( |B3D requires the computation of the following integral : 



J, 



nl Cl, . . . , tn) 



/+00 r r+OD r+oo 

dy P« / ... / ds2... dsn F{x, y,0,S2,... , s„) 
oo L Jo Jo 

f+oo r+oo 
+ Vl^Jg ... / dsidS3...dSn F{x,y,Si,0,S3... ,Sn) + ■■■ 

Jo Jo 

/•+00 f+oo 
+ ^reg / •••/ dsi . . . dSn-l F{x,y, Si, . . . , Sn-1,0) . 

Jo Jo 



(B9) 



This expression corresponds to a sum of n trajectories ; the j trajectory having one Keller bounce 
on apex j (with a regular part of the diffraction coefficient noted P^eg, j = 1 . . . n) and Kirchhoff 
contributions from the other apexes (for instance Eq. (|T4|) corresponds after transverse integration to 

J2ie,i)). 

We obtain 



J2(4,4) = (/?^i), + /?^?J 



„ikL+3in/4 



'"' '''' AkVs^M: 



(BIO) 



and 
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^3(^1,^2,^3) 



^ikL+3in/4 



SkVSTTkL 



T^ill 1 + - arctE 



reg 



IT 



^3^ 



+ vi% h + ~ ^^^*^ 



his 



+ Vill 1 + - arctE 



reg 



IT 



The expression of J„ when all the £'s are equal is : 



•^nV'^j •••, ^-j 



pikL—inir/i 



^ ^req •Jn—1 



where 



«^n 



+ 00 

dXi...dXn en^?+(^2-xi)2+. 



. + (x„-x„_i)2+z2 



^2^ 



(Bll) 



(B12) 



(B13) 



It is shown in [Appendix C| that J„ = (e*'^/^A/vr)" (n + 1) ^/^. From this result and Eq. ( |B12| ), 
formula (|l5|) follows immediately. 



• In section ^ in order to compute the first order correction to the contribution of a family whose 
boundary partly coincides with the frontier of the billiard, one needs to compute the following integral : 



M2(h,i2) = Wl,i2) - {2ikfD2 



dy I dsi / ds2 exp ^ — 
-oo Jo Jo I ^ 



si -yf , {s2 + sif , {y-s2 



+ 



+ 



'ii- X 
(B14) 



This equation corresponds to the transverse integration of (pH) after the removing of the contri- 
bution of the the direct path (i.e. of Go). The last integral in the r.h.s. of ( B14| ) corresponds to the 
orbit going from r to r and bouncing on the boundary of the family which is also a frontier of the 
billiard. The term 12(^1,^2) corresponds to the direct diffractive trajectory. One obtains 



M' 



2\.i^l,i^2) 



pJ,kL+in/2 

AnkL 



y^£^ + L avctg J j- 



(B15) 



Here we want to develop a point stated in the main text : if f lies near the part of the optical 
boundary that coincides with the frontier of the billiard, the Green function has four contributions 



which, after integration over y, give the same contribution as ( B15 ) 
If one defines 
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h^{x,y,si,S2) = exp<{y 
/i^(x,y,si,S2) = exp<J — 



h'^(x,y,si,S2) = exp<J — 



{s2-y)'^ {si-S2y {y-si)"^ 

X £i (-2 — X 

{s2 + y? ^ {si-s2f ^ jy-sif 

X £i (-2 — X 

{s2 + y? , {si-s2f , {y + sif 



+ 



+ 



X 



t2-X 



h'^{x,y,si,S2) = exp<J — 



{s2-y? ^ [81-32^ ^ jy + si)^ 

X t\ ^2 — X 



(B16) 



then, for a point r near the boundary of the orbit which coincides with the frontier of the bilhard, the 
four contributions to the Green function integrated transversely to the direction of the orbit read : 



/•+00 r ^+oo r+oo p+oa p+oo 

M2{x) = {2ikfD2{x,£2,£i) dy / dsi ds2h^{x,y,si,S2) - dsi ds2h^{x,y, 81,82) 

Jo Uo Jo J-00 J-00 



(B17) 



and 



f+00 f+oc f+00 

M^ix) = i-iy+^i2ikfD2ix,i2,h) dy dsi ds2h^ix,y, 81,82) 

Jo Jo Jo 



for j = 2,3,4. (B18) 



Note that the transverse integration (over the variable y) is only possible here from to +00 
because one is near the frontier of the billiard (see Fig. |2^ ). The four contributions (B17, BT8| ) 
correspond to different paths going from r to r : M2 corresponds to a path going from r to 82, to si 
and back to r (for this part one has to withdraw the semi-classical Green function) ; M2 corresponds 
to the path going from r to 5*2 with a reflection on the boundary of the orbit which coincides with the 
frontier of the billiard, then going from 82 to si and back to r etc ... This is illustrated in Figure p3. 



The Mg's separately depend on x, but one obtains 



Ml (x) + M^ [x) = ^ , ^ VfT^ and M^ (x) + M| (x) = — j-^— arctg \ / -^ 



AnkL 



A-nk 



(B19) 



Hence, when r lies near the frontier of the billiard, one obtains M2 (x)-|-M2 (x)-|-M2 (x)-|-M2 (x) = 
M2, where M2 is given by ( B15| ) and corresponds to transverse integration when r is not in vicinity 
of the frontier of the billiard. Thus the result of the transverse integration of the Green function does 
not depend on the position x of the point r along the orbit. This directly leads to formula (|2l|). 



• In Section ^, for evaluating the contribution of a diffractive orbit jumping from one boundary of 
a family to another one needs to compute the following integral (which is the transverse integration 
of Eq. (H)) : 
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(3) 




(1) 




(4) 



Figure 24: Schematic representation of the path encompassed in the contribution M^ (Eqs. ( B1? , B18| )) to the 
transverse integration of the Green function (j — 1...4). The plot labelled (j) corresponds to M^ix). One has 
here four different contributions because the initial point r lies along the part of the boundary of the family 
that coincides with the frontier of the billiard. The simpler case that r is not in vicinity of the frontier of the 
billiard is represented in Fig. nG. 



/+CX) K r+oo 

dy / dsi / ds2F{x,y,si,S2). 
oo J —oo Jo 

The y integration is trivial and (as now usual) removes the x dependence. One obtains : 
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y/SnkL 

where L^ is the length of the difFractive orbit close to the family of length L and A = L^ — L 
(A ~ C,"^ L / {2iii2)) ■ In the above expression, the last term simply re-expresses the previous one 
using the modified Fresnel integral defined in Eq. ( [A2[ ). The longitudinal integration (which simply 



amounts to a multiplication by L) of ( B21 ) yields the contribution to p{E) of the configuration studied 
in Section y. More precisely, this contribution reads — (L/vr) Im A'^2(^i)^2)- This results in Eq. (pq). 



35 



Appendix C 

The purpose of this Appendix is the exphcit computation of the integral ( [B13| ) 



cxD r+oo 



/O ^0 

where ^n{x) is the following quadratic form 



dxi...dxn e 



i^n{x) 



(CI) 



$„(x) = x\ + (Xi - X2f' + ... + {Xn-l - Xnf + X^ • 



(C2) 



Note that in all this Appendix we denote n-dimensional vectors x = {xi, ...,Xn)- A key point 
in the evaluation of integral ( |C1| ) is the existence of a group generated by a set of transformations 
\Tj}^^ .^^ which leaves the quadratic form invariant : 



Tj{x)=^ with 






(C3) 



where j = l...n and we have adopted the convention x_i = Xn+i = 0. 

These transformations are inversions (T? = 1) and they generate a finite group (of the An type, 
see e.g. |Q). We give below a method of calculation of ( pl[ ) which does not require knowledge of the 
theory of finite groups. 



The quadratic form 



can be naturally rewritten in the form : 



$„(x) = 2 ^MijXiXj 



i=l 



where the n x n matrix M is the following : 



(C4) 



M 



V 



1 


-1/2 

















1/2 


1 


-1/2 

















-1/2 


1 


-1/2 .. 























-1/2 


1 


-1/2 

















-1/2 


1 



/ 



(C5) 



Due to the simple tri-diagonal structure of M_, it easy to perform a Gaussian decomposition of the 
quadratic form <!>„ recursively. This results in 



re-l 



fe=l 



^'kiXk - TTl—) + A„X„ . 
ZAk 



(C6) 
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where 



Afe = ^ k = l...n . (C7) 



From (C6,P7|) it is clear that the determinant of M_ is 



71-1-1 

detM = Ai...An = ^. (C8) 



Let us now introduce n vectors V^ such that 



Mij = V\V^ . (C9) 

From ( |C5| ) it is clear that (if these vectors exist) they are of unit length and that the angle between 
different vectors equals either 7r/2 or 27r/3 (the cosines being either or —1/2). One possible solution 
of Eqs. ( |C9|) can be written in the following form : 



v} -- 


= (yAT,0,...,0) , 


y2 -_ 


= (^y^^^VA^^O,.. 


v^ = 


^ ^^'-'SyA^ 


V"^ = 


= (0'-'0':r^ 



2v An_i 



Afc,0,...,0) 
Aji) 



(CIO) 



where the A^'s are defined in ( iCqJCTl) . From geometrical considerations it is clear that all other 
solutions can be obtained from Eq. ( |C10[ ) by applying overall rotations (and inversions). 
Using the vectors V] it is possible to define new coordinates (y^) from the relation 



]L = Y.'^,VK (Cll) 

Note that 



M = J^ ^^ ^jV\V^ = $„(x)/2 . (C12) 

The Jacobian of the transformation (|C11| ) is 
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(C13) 



It can be computed from Eqs. ( Cl[l| ) or by taking the square of expression ( C13 ). This yields 
J'^ = detM. Therefore 



oo r+oo 



1 



|detM|i/2 



... / dyi...dyn e 



2i{yl + ... + yl) 



(C14) 



where the integration is taken over the interior Vt of the hyper-spherical simplex defined by the n 
vectors V^ , 



n 



/. ^3 y^ ^^^ a^j > ; J = 1, ..., n 



(C15) 



i=i 



• The next and final step of the computation of ( |C1| ) is to show that the integration domain Q, in 
( |C14 ) is a relatively simple sub-part of the whole n-dimensional space. To understand geometrically 
the structure of this region it is convenient to add a new vector y""*"^ to the list (P1C|) such that 



V^ + ... + V'' + V 



n+1 



. 



(C16) 



Angles formed by V^~^ with the other W's are straightforwardly obtained from the expression 



1/2, y"+^p = for j = 2, ..., n - 1 and ]/"+^y"+i = ^.^^. y\p 



(|C5| ) : V^+^.V} =Y^^^.V^ -- 
1. 

From the (n + 1) vectors V}, ...jF""*"^ one can define (n + 1) regions Qj : 



n^ 



V^ Xj W^ and Xj > ; j = 1, ..., 



n 



i=i 



(C17) 



where the n vectors Wj include all (n + 1) vectors V^ but the vector V^ . In these notations the region 
O in ( |C15| ) coincides with ^n+i- 

It is almost evident that these (n + 1) regions cover the whole y space without common intersection 
points. A formal proof of this statement can be the following : 

An arbitrary point y of the n-dimensional space has a unique decomposition on the non-orthogonal 



basis of the V^-^'s {j = 1, ...n) as given by ( Cll ). If all Xj > then y G i^n+i, otherwise the set of 
coordinates Xj is divided into 2 sets x'^ and x'L of positive (x^ > 0) and negative {x'L < 0) coordinates. 
Denoting by z^ = —x'^ and by z^ the maximum of the z^'s we get : 






(C18) 
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Using the definition ( C16| ) one can express VJ as a function of the other V_^s and rewrite the above 
expression as 

y = 5^(x„ + z^)Z" + E(^7 - ^p)y^ + ^lYT"-' ■ (C19) 

As z^ > zp all coefficients in this sum are non- negative and the point y hence belongs to Q^. It 
is thus clear that regions Jlj's have no common points except on boundaries where some Xi = 0. It is 
also clear that the union of all the rjj's {j = l...n + 1) covers all space since, given an arbitrary point 
y, one can assert unambiguously from the above procedure at which of the Oj's it belongs. 

• Each region Qj is defined by n vectors W^ obtained from the (n + 1) vectors y_ s by ignoring 
Vj' . The convenient rearrangement of vectors W^ 's is the following : 

(li:^)i>,>„ = {V'^\V^'-',.:,Yr,Vr^\v\V^...,V^-') . (C20) 

Let us now construct the matrix of mutual projections Nij = W^ .WK It is easy to check that this 



matrix coincides with the matrix M defined in Eq. (C5). Therefore the vectors W^ {j = l...n) will 
have the same mutual positions as our initial vectors V^-^'s. As we noted above, it means that region 
Qj for all j {j = l...n + 1) can be obtained from the initial region 17 (= r2„+i) by overall n-dimensional 
rotations (and possibly by inversions). But the integrand in the r.h.s. of Eq. ( pl4| ) is invariant under 
such transformations, therefore its integration over any of the J7j's is the same and : 






n n + 1 J_^ 7_oo n + 1 



Combining this result with Eqs. ( C14| ) and (C8) one obtains 



Jn = n.. rix,...dxn e^ ^«(^) = ^-f^^ . (C22) 

Jo Jo [n + ^r'^ 

which is the result needed for the explicit computation of expression ( pi2| ). From the approach 
exposed in this Appendix, one can state the more general result : 

/ ... / dxi...dxnf{^nix)) = -—-^j^—— / drr^-'f{r'). (C23) 

Jo Jo {n + iy/^T{n/2) Jq 

Appendix D 

In this Appendix we derive the explicit form of I„(^, ...,£) defined in ( [B3D , or equivalently we give 
the value of the coefficients Cn appearing in Eq. (p8|). These computations extend the results of 
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[Appendix B| (Eqs. ( |B4| ) to (|B7])) and are valid for any n. However, they are restricted to the case 

4 = ... = 4=^. 

For evaluating the integral ( |B3D it is customary to make several manipulations : one performs 
the y integration in the first term of the r.h.s. In the second term, one can decrease by one the 
number of variables of integration easily, since this term is simply an elaborate manner of writing 
G(y) Go{r,r,E). After a scaling on the variables {yi = Siy^k/{2£)) one obtains a result which can be 
cast in the form : 



where 
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(Dl) 



(D2) 



and ^(x,y) is the following quadratic form (we denote n-dimensional vectors y = (yi, ...,y„)) 



^(x, y) = {x- yif + (yi - 2/2)^- + {Vn-i - Vn? + {Vn - xf . (D3) 

/„ in ( p2| ) can be expressed simply in term of the function ipnix) defined by : ' 

/> + 00 I' + OO 

Mx)= / ... / (iyi...(iy„e^*(^'M). (D4) 

Jo Jo 

One first notices that for x large and positive, one can neglect the boundary effects in the integral 
( p4D defining ipn{x) and thus : 



1pn{x) > Vn(+00) 

x^+00 

Hence (D2) can be written as : 
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+00 p+oo 

In= I dx [ll^ni+Oo) - ll^nix)] = / dx X 1p'n{x) . 

Jo 



The function ip'^ in (D6) can be cast in the form 



(D6) 



" p+oo p 



+00 p+oo 

dy,...dyn6{y^)e'^^'''yJ 



(D7) 



This is done by first changing variables in ( p4| ) (y^ = x -\-tj), then deriving with respect to x and 
finally coming back to the original variables yj. 



40 



Inserting this expression in ( |D6| ) and renumbering the variables in the integral, one obtain the 
following expression for /„ : 



in = ^ {Vm) where (y„ 



m,=l 



+00 r+00 



^0 



dyi...dyn Vm e 



i^n{y) 



(D8) 



and the quadratic form <!>„ is defined in Eq. (|C^). The n integrals {ym) are computed by means of 
the auxiliary integral Pm defined as : 



i^n{y) 



00 /■+OO Q^ 

... / dyi...dy, 
JO *^ym 



00 r+00 



dym JO JO 

-Ai{ym) + 2i{ym+i) +2i{ym-i) • 



dyi...dyrn--idym+i-dyne " - '!/™=0(D9) 



(DIO) 



with the convention that (y-i) = {yn+i) = 0. The integrals P^ are easily calculated using the 

results of [Appendix C and noticing that <l>„(y) _ is a sum of two independent quadratic forms 

— — 'ym — 

$m-i(yi) ■■■lUm-i) and (^n-m{ym+i-, ■■■,yn)- Therefore the integral of the r.h.s. of ( |D9| ) is the product 
of two integrals of type (|Cl] ) (whose explicit form is given in ( C22|) ) : 



-* m — •Jm—1 'Jn- 



(e^-/^^) 



n-l 



(n — m + 1)]'^/^ 
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(Dll) 



Then it is a simple matter to solve recursively the system of equations formed by ( D10| ) and to 
express the (2/m)'s in term of the Pm's. This yields : 



{y-n 



^ 9=1 



2(n + l)^ 
^ ' 5=1 



(n + 1 - g) Pg 



m - q) Pq 



and 



(D12) 



in = Yj (y- 



m=l 



gi7r/2 '^ 

m/^^Y.^qin-q + l)]Pq. 

9=1 



Using Eq. ( Dll| ) one obtains the final formula 



(D13) 



pin/2 ■"■ 



1 \/(l{n-q + l) 



(D14) 



From Eqs. (B8,D1,D14), the coefficient C„ appearing in Eq. ([171) reads 



41 



n—1 



1 1 



which coincides with the results obtained for n = 2, 3, 4, 5 in Appendix B| by a different method. When 
n — > oo the sum over q can be substituted by an integral and : 

1 r^ dr 1 

hmC„ = -L/ -r=^ = ^- (D16) 
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